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ABSTRACT 

Analysis is made of the stresses around a rigid circular inclusion with cemented 
boundary in a rubber sheet which is under simple tension at a distance from the 
inclusion. The analysis is performed by the photoelastic method, in which data 
are obtained from observations of the changes produced in a beam of polarized 
light passing through the transparent model. The fundamental principles of the 
method and a detailed account of their application are given. The possibilities of 
the use of rubber as a model material are demonstrated. The stress-optical coeffi- 
cient of a certain soft vulcanized rubber compound is found to be 2,030 brewsters. 
This is in agreement with the results of previous work. The maximum stresses 
near the boundary of the disk are found to be about 50 percent greater than the 
average applied stress. The experimental results are compared with those given 
by a theoretical treatment of this problem, and quite satisfactory agreement is 
obtained. 
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I. INTRODUCTION 

The problem of the stress distribution around a rigid inclusion in a 
body under tension is of considerable interest in a number of different 
fields, and yet there has been comparatively little work done on it. 
The related problem of the stresses around a hole has received much 
more attention. In the case of rubber, the problem of the rigid 
inclusion is particularly important not only with respect to inclusions 
of appreciable size, but also with respect to the microscopic and sub- 
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microscopic filler particles, which must give rise to some rather com- 
plex stress distributions. The present investigation of the stresses 
around a circular inclusion in rubber was originally undertaken with 
the aim of solving the much more complicated problem of the distribu- 
tion of stress in rubber containing filler particles. Some preliminary 
work was done, on the mutual effect of the stresses arising from the 
presence of more than one inclusion, but the investigation was not 
completed. 

For the experimental solution of such stress problems the extremely 
useful method of photoelastic analysis has come into prominence in 
recent years, especially since the publication by Coker and Filon of a 
comprehensive treatise l on the subject. The necessary data are 
obtained by observing in polarized light the patterns obtained when a 
transparent model is stressed. 

Many years ago Maxwell 2 outlined the theory underlying the general 
methods of such analysis, but it remained for Mesnager, 3 Coker, 4 and 
others to develop practical procedure 5 6 and the materials 7 for the 
models. Extensive use of the method has been made in the study of 
engineering materials with models constructed of glass, Celluloid, 
Bakelite, and other substances, especially the organic plastics. 

Since the results of a previous investigation 8 of the photoelastic 
properties of soft vulcanized rubber have indicated that some rubber 
compounds fulfill many of the requirements for a desirable model 
material, it was decided to study the problem of the stresses around 
a rigid circular inclusion with a cemented boundary in a rectangular 
model consisting of a sheet of transparent, soft vulcanized rubber 
under tension in one direction. In this manner, therefore, contrary 
to the usual practice, the problem could be studied with a model of 
the very material for which the solution was desired. Furthermore, 
by a definite illustration, it was desired to give a demonstration of 
the possibilities of the use of rubber as a model material in photo- 
elastic work, and to present in convenient form for those unfamiliar 
with this field a detailed outline and illustration of the exact experi- 
mental procedure used for the solution of problems by this method. 

II. THEORY OF STRESSES 
1. PRINCIPAL STRESSES 

Through every point in a plane section under stress there is a pair 
of mutually perpendicular lines across which the stresses are normal, 
the shearing stresses being zero in this case alone. These normal 
stresses are known as the principal stresses at this point. Across all 
other lines through the point, in general, shearing stresses exist. The 
directions of the principal stresses, if joined in continuous curves 
throughout the body, form a system of mutually orthogonal lines. 

i E. G. Coker and L. N. G. Filon. A Treatise on Photoelasticity (Cambridge Univ. Press, 1931). 

>J. C. Maxwell. On the equilibrium of elastic solids. Trans. Roy. Soc. Edinburgh 20, 87 (1853); Re- 
printed in Scientific Papers 1, 30 (Cambridge Univ. Press, 1890). 

8 A. Mesnager. Contribution to the study of the elastic deformation of solids. Ann. Ponts et Chaussees 4, 
128 (1901); Determination of stresses in a glass model. 16, 133 (1913). 

« E. G. Coker. Photoelasticity for engineers. Gen. Elec. Rev. 23, 870 (1920); 23, 966 (1920); 24, 82 (1921); 
24, 222 (1921); 24, 455 (1921). 

* A. L. Kimball. Optical determination of stress in a transparent material. J. Opt. Soc. Am. 5, 279 (1921) 
«H. B. Maris. Photoelastic investigations. J. Opt. Soc. Am. 15, 203 (1927). 

1 R. B. Carleton. Suitability of materials for photoelastic investigations. Rev. Sci. Instr. 5, 30 (1934). 

• Wilfred E. Thibodeau and Archibald T. McPherson. Photoelastic properties of soft, vulcanized rubber 
J. Research NBS 13, 887 (1934) RP751. Reprinted in Rubber Cheni. Tech. 8, 183 (1935). 
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Knowledge of the magnitude and direction of the principal stresses at 
every point in the body furnishes a complete solution of the stress 
problem, both the normal and shearing stresses in any other direc- 
tion being derivable by means of simple equations. 

The symbols P and Q will be used to represent the maximum and 
minimum principal stresses, respectively, at a point in the plane of 
stress. Since neither normal nor shearing stresses can exist on a free 
boundary, the value of one of the principal stresses at a free boundary 
is zero. Tensile stresses will be represented by positive values for P 
or Q, while compressive stresses will be represented by negative values. 

2. ISOCLINIC LINES 

As an aid to the expression of the directions of the principal stresses 
at all points in a plane section under stress, it is convenient to define 
a family of curves, called isoclinic lines, along each of which the di- 
rection of the principal 
stress P makes a con- 
stant angle with the r^\. ^ 

direction of applied ° ^^- _ 0*A0 j _^- — — 0+A0 

stress. This constant 
angle, the parameter 
of the isoclinic, here 
denoted by <£, is the 
same as that between 
the normal to the 
direction of the ap- 
plied stress and the 
direction of the Q 
stress. The family of 

, 4 ' Figure 1. — Method of plotting isostatic lines from the 

serves to specify com- isoclinic lines. 

pletely the directions 

of the principal stresses at all points in the plane section. It will be 
shown later that the isoclinics in a transparent medium under stress 
may become visible in polarized light under the proper conditions, and 
can be traced from direct observation. 

3. ISOSTATIC LINES 

Another family of curves, called isostatic lines or simply lines of 
stress, can be drawn so as to indicate by their direction at any point, 
the direction of one of the principal stresses. In other words, for 
example, the direction of the Q stress at any point is given by the 
tangent to the Q isostatic at that point. 

This manner of plotting the isostatics from the isoclinics is shown 
in figure 1. Let <j> and <£+A<£ be adjacent isoclinics. </> then is the 
locus of all points at which the principal stress line makes an angle <j> 
with the zero position, and at the points along <j>-\- A<£ the principal 
stresses are inclined at the angle 0+A0. If now at certain intervals 
along <£, lines are drawn making the angle $ with the zero position, 
these will be tangent to the principal stress lines. These lines are 
continued to a point midway between the two isoclinics, and from there 
are continued through 0+A<£ at an angle </>+A$. A continuation of 

44637—38 10 
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this process will furnish an envelope to which the principal stress lines 
can be drawn. After one set of principal stress lines has been com- 
pleted in this manner, the other family of stress lines can be drawn 
in at right angles to the first or it can be obtained from the isoclinics 
in the same manner as the first. 

The isoclinics and one set of isostatics, the Q lines, for the present 
case are shown in figure 2. These lines in themselves are not suffi- 
cient for a complete solution of the stress problem, since this requires 
the magnitudes of P and Q as well as their directions. 




APPUED 
TENS/ON 



/SOCL/N/CS 

fSOSTAT/CS 

Figure 2. — Isoclinic and Q-isostalic lines. 

4. MAXWELL'S EQUATION OF EQUILIBRIUM 

Maxwell, in his fundamental work, presented a pair of simplified 
equations of equilibrium: 9 



P-Q 


dP_ 


P-Q 


h ds P 
i>Q_ 


r r 


ds„ 







(la) 
(lb) 



On integration eq la yields: 



P=P, 



H 



S >P-Q 



ds r 



(2) 



• J. C. Maxwell. On the equilibrium of elastic solids. Trans. Roy. Soc. Edinburgh 20, 87 (1853); Reprinted 
in Scientific Papers I, 30 (Cambridge Univ. Press, 1890). For a derivation of these equations, see page 129 
of S. Timoshenko, Theory of Elasticity (McGraw-Hill Book Co., New York, N. Y., 1934). 
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in which P and Q are the principal stresses, s P is the distance along 
the P line from the point where P=P , S T is the value of s P at the 
point considered, and r Q is the radius of curvature of the Q line at 
this point. Similarly, integration of eq lb gives: 

rs*p_n 

Q=Qo+j V^ s « (3) 

in which s 9 is the distance along the Q line from the point where 
Q=Qo, S Q is the value of s Q at the point considered, and r P is the 
radius of curvature of the P line at this point. P and O > the con- 
stants of integration, must be separately determined. 

It is to be noted that these equations are merely equations of 
equilibrium, and so involve no elastic constants. Hence they will 
apply as well to plastic as to elastic substances ; and this includes sub- 
stances which have been stretched beyond the elastic limit and mate- 
rials like rubber, which do not conform to Hooke's law. The only 
restrictions are that the material be homogeneous and in a state of 
equilibrium. 

In practical work with a plot of isoclinics and isostatics it is some- 
what difficult to measure directly the values of r Q and r P , and it is 
easier to compute them from the relation r= arc/angle =AaS7A$, where 
AS is the arc on the stress line normal to the path of integration and 
bounded by the isoclinics of parameters <f> and </>+A0 in the neigh- 
borhood of the point. 

Furthermore, there are a number of advantages of improved 
applicability in expressing results in terms of ratios of the principal 
stresses to some applied stress rather than in terms of the principal 
stresses themselves. In the first place, the results are then inde- 
pendent of the stress units employed. Secondly, the results are 
immediately applicable to cases in which the applied load is different 
from that which happened to be employed in the particular case 
under consideration. Finally, a comparison with the results of a 
theoretical treatment is possible at each point of the sheet in terms of 
merely the coordinates of the point. 

Therefore, eq 2 and 3 are modified by a substitution for r and by 
dividing both sides by the applied stress T, defined in the present case 
as the ratio of the applied load to the reduced cross-sectional area of 
the model. 

T T Jo T Asq F w 

<3_&, c**p-qh , ,,, 

T~ T^ Jo T As P aSa w 

An alternative modification of eq 2 and 3, suggested by Filon 10 
gives: 

P P C* p Q 

T = T° + I ~~T cot ^* ( 6 ) 

r = T~I!^r^ cot ^ (7) 



i° L. N. Q. Filon. Graphical determination of stress from photoelastic observations. Engineering 116, 511 
(1923); Reprinted in Brit. Assn. Rep. 91, 350 (1923). 
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Here \I/ P and fa represent the angles through which the P line of stress 
and the Q line of stress, respectively, must be rotated in a counter- 
clockwise direction to bring them upon the isoclinic, and <j> is the value 
of 4>, where P=P and Q=Q . 

The photoelastic method, as shown below, makes possible the 
evaluation of the quantities either in eq 4 and 5 or in eq 6 and 7. The 
question of which method to use is a matter of judgment and depends 
on the configuration of isoclinics and isostatics at the point considered. 
The ideal case for eq 4 and 5 is where ^=0, that is, where the line of 
stress is itself an isoclinic. The ideal case for eq 6 and 7 is that in 
which the isoclinics are inclined at an angle of about 45° to the lines 
of stress and are closely packed. If \f/ is too small the errors will be 
large in this method. In the present work both methods were used, 
sometimes even along the same line of integration. 

III. PHOTOELASTIC THEORY 
1. STRESS-OPTICAL COEFFICIENT 

When light passes at normal incidence through a sheet of isotropic 
transparent material stressed in its plane, it emerges, as in the case of 
doubly-refracting uniaxial crystals, as two plane-polarized rays called 
ordinary and extraordinary rays, one polarized in the plane of one of 
the principal stresses and the other polarized at right angles to it, in 
the plane of the other principal stress. The velocities of these two 
rays are different, and depend on the intermolecular and interatomic 
forces in the directions of the principal stresses. The velocity of the 
ordinary ray is found to be proportional to the P stress plus a con- 
stant, and the velocity of the extraordinary ray to vary in exactly the 
same way with Q. The difference of the velocities is then proportional 
to the difference of principal stresses. The difference in the velocities 
of the two rays introduces a phase difference as they emerge from the 
medium. The phase change in each ray, as it passes through the medi- 
um, is proportional to the thickness of the medium and inversely pro- 
portional to the velocity of that ray. The phase difference, or retarda- 
tion of one with respect to the other, is then proportional to the 
thickness and to the difference of the reciprocals of the velocities. 
Since the difference in the two velocities is extremely small compared 
with the velocities themselves, the difference of the reciprocals is 
approximately proportional to the difference of the velocities them- 
selves. 

It can thus be seen that the retardation of the ordinary with 
respect to the extraordinary ray at any point in the body is propor- 
tional to the difference of the principal stresses P and Q and to the 
thickness of the sample. That is, 

R=C (P-Q) t, (8) 

where R is the relative retardation, t is the thickness of the medium, 
and (7, thus defined, is a constant known as the stress-optical coeffi- 
cient. 11 When the relative retardation is expressed in angstrom units, 
the stress in bars (1 bar=10 6 dynes/cm 2 =1.0197 kg/cm 2 = 14.504 
lb/in. 2 ), and the thickness in millimeters, the stress-optical coefficient 
is given in brewsters. A brewster is thus 10~ 13 cm 2 /dyne. 

"E.G. Coker and L. N. G. Filon. A Treatise on Photoelasticity, p. 185 (Cambridge University Press, 
1931). 
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When the value of C for the material is known, the value of (P—Q) 
at any point in the medium can be obtained by a determination of 
the retardation and a measurement of the thickness. This law, as 
expressed in equation 8, has been fairly definitely established for 
materials like Bakelite, Celluloid, and glass, and more recently for 
some vulcanized rubber compounds. 12 A compound (compound B, 
composition given below) which conformed to this relation was selected 
for the present investigation. 

2. LOCATION OF ISOCLINIC LINES 

It will now be shown how, as previously stated, in polarized light 
under the proper conditions the isoclinic lines become visible and can 
be traced from direct observation. Light which has been polarized 
by passage through a nicol prism, called the polarizer, passes through 
the sheet and is observed through a second nicol prism, the analyzer. 
At all points where the vibration plane of the polarizing nicol is at right 
angles to the direction of one of the principal stresses, one of the usual 
two plane-polarized rays passing through the medium is absent, since 
there is then no component of vibration corresponding to the plane of 
vibration of this ray. Under the same conditions if the analyzing 
nicol is set at right angles to the polarizing nicol, the second plane- 
polarized ray, which is actually present in the specimen, is not observed 
since it can not pass the second nicol. That is, as far as the point in 
question is concerned, the result is the same as if the model were not 
interposed. Hence a dark band or brush will appear marking the 
locus of all points in the model where the principal stress directions 
correspond with the vibration directions of the nicol prisms. In order 
to obtain a complete map of the isoclinics it is necessary merely to 
rotate the nicols by stages through 90°, keeping their vibration planes 
at right angles and tracing out one resulting isoclinic at each setting. 
The orientation of the prisms at each position is the parameter of the 
resulting isoclinic. 

IV. EXPERIMENTAL ARRANGEMENT 
1. APPARATUS 

The apparatus used for the photoelastic determinations consisted 
essentially of a light source, two nicol prisms, means for supporting 
and stressing the model, suitable lenses, and a screen. The light 
source was a 1,000-watt projection lamp, with a water filter for the 
reduction of heating effects. By means of condensing lenses the light 
rays were brought to a focus within the nicol prisms and were reduced 
to parallelism in the region between the nicols where the model was 
placed. A further refinement, the placing of divergent lenses on either 
side of the nicol prisms so as to transmit the rays through them in a 
parallel beam, was found to be unnecessary with the large nicol prisms 
used. Both prisms were equipped with divided circular scales so that 
their inclinations could be recorded. 

An arrangement of two displaceable quartz wedges known as the 
Babinet compensator was inserted between the model and the analyz- 
ing nicol when it was desired to measure the retardation produced in 

i» Wilfred E. Thibodeau and Archibald T. McPherson. Photoelastic properties of soft, vulcanized rubber. 
J. Research NBS 13, 887 (1934) RP751. Reprinted in Rubber Chem. Tech. 8, 183 (1935). 
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the stressed model. The compensator itself produced a retardation 
in the opposite direction amounting to 1 1 ,400 angstrom units for each 
centimeter of displacement of the wedges. In this manner the net 
retardation could be reduced to zero, and measurement made of the 
amount of retardation required to accomplish this. The compensator 
was mounted so that it could be turned until its axes were parallel to 
the principal stresses at the point under observation. The final lens 
focused an image of the model on a piece of plain white paper, which 
served adequately as a screen. 

2. MODEL 

The model was a rectangular sheet of soft transparent rubber 46.0 
by 29.8 by 0.651 cm, in the center of which was placed a disk of stiff 
rubber-flooring material, the moduli of which were very great in 
comparison with those of the soft rubber. The disk was 5.22 cm in 
diameter and 0.62 cm in thickness. The periphery of the disk was 
coated with rubber cement (Vulcalock), and after vulcanization the 
disk became firmly embedded in the model. 

The composition of the soft rubber was as follows: 

Parts by weight 

Rubber 100. 

Sulphur 2. 

Zinc oxide (Kadox brand) 1. 

Stearic acid 1. 

Tetramethylthiuram disulphide 0. 3 

This was cured for 45 minutes at 125° C in a steam press. The 
weight of the model was 858 g. It was loaded vertically by means 
of wooden clamps bolted to each end of the model. The weight of 
the lower clamps and bolts was 170 g. 

V. PROCEDURE 
1. DETERMINATION OF STRESS-OPTICAL COEFFICIENT 

For the purpose of determining the photoelastic characteristics of 
the rubber, a separate tensile specimen of flat dumbbell shape was 
prepared from the same compound, and stretched. The results of 
observations in the central region of uniform stress are shown in 
table 1. Since here Q,=0 and P=T, by the application of eq 8 to 
this case, one finds: 

R=CTt, (9) 

where T is the applied stress, measured on the actual cross section 
and t the actual thickness. If a is the compensator constant, the 
retardation per unit displacement of the compensator, and X the 
displacement, 

R=aX=CTt 

T-lv (10) 

1 ~Ct X 
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Table 1. — Data for calculation of stress-optical coefficient 

[Original thickness 0.278 cm; original width 0.632 cm] 

W= Applied load. 

X= Distance between gage marks. 

X= Displacement of compensator. 



w 


L 


X 


kg 


cm 


cm 


0.000 


4.98 


0.000 


.055 


5.13 


.145 


.123 


5.31 


.330 


.173 


5.46 


.485 


.223 


5.62 


.650 


.273 


5.77 


.815 


.323 


5.98 


.970 


.373 


6.18 


1.140 


.423 


6.40 


1. 345 


.473 


6.63 


1.510 



For pure gum rubber vulcanizate, which was used here, Poisson's 
ratio has been found 13 M to be 0.500 for elongations in the range under 
consideration. For an isotropic medium with a Poisson's ratio of 
0.5, it may easily be shown that 

r-r,(i+.) (ii) 

and 

"7TT< (12) 

where T is the applied stress calculated on the original cross section, 
t the original thickness, and e the strain, or fractional elongation. 
Then 



m+ t )- 



w 



A r Vl+« 



(13) 



This equation predic ts pro portionality between the quantity T (l+€) 
and the quantity X-yjl-^-t. 

A graph of the experimental data of table 1 verifies this propor- 
tionality and has a slope 15 of 2.06 kg cm" 3 . On multiplication of this 
slope by the thickness, the value of v/Cis found to be 0.572 kg cm" 2 . 
It is not necessary to compute the actual value of C. 

Turning now to a consideration of the stretched model with the 
circular insertion, the mean thickness t m of the model being 0.611 cm, 
one can apply eq 8 to obtain 



P-Q- 



o~ 
"Ctr 



X= 



0.572 
0.611 



^-0.936A r 



(14) 



Since the thickness of the model varied from point to point, strictly 
the individual values should be used instead of a single mean thick- 
ness. This more rigorous method was applied at a number of selected 

M Wilfred E. Thibodeau and Archibald T. McPherson. Photoelastic properties of soft, vulcanized rubber 
J. Research NBS 13, 887 (1934) RP751. Reprinted in Rubber Chem. Tech. 8, 183 (1935). 

i* William L. Holt and Archibald T. McPherson. Change of volume of rubber on stretching: Effects of time, 
elongation, and temperature. J. Research NBS 17, 657 (1936) RP936. Reprinted in Rubber Chem. Tech. 
10, 412 (1937). 

16 From this slope a value of the stress-optical coefficient C may be calculated. Here to =2. 78 mm, and 
<r= 11,400 A/cm. 

C ° (2.78) (2I06H98O) 1Q-3 " 2 ' 030 brewsters. 

This result is in reasonable agreement with the value 2,096 brewsters, the average obtained over a much 
larger range of stresses for seven samples of this compound in the previous investigation by Thibodeau 
and McPherson (reference 13). 
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points, and was found to yield results differing from the simpler method 
by only about 5 percent at the most. In view of the difficulties of 
measurement and the increased complexity of calculation, it did not 
seem worth-while to utilize this more refined procedure. 

The load applied near the center of the model included, in addition 
to an external load of 26.75 kg, the weight of the lower clamps and 
half the weight of the model. When this total load of 27.35 kg was 
applied, the distance between the clamps increased from 41.2 to 46.1 
cm, the width of the model at the center decreased from 29.8 to 28.3 
cm, and its mean thickness decreased from 0.651 to 0.611 cm. Dividing 
27.35 kg by 17.3 cm 2 , the reduced cross-sectional area, one obtains 
1.58 kg/cm 2 as the value of T, the applied stress. Finally, 

P^_0,936^ a593Z (lg) 



T 1.58' 

This is the relation employed in the calculations involving the data 
obtained with the compensator. 

2. LOCATION OF ISOCLINIC LINES 

With the plane of vibration of the polarizing nicol vertical and the 
plane of vibration of the analyzing nicol horizontal, the external load 
of 26.75 kg was applied to the model along the vertical direction. 
The isoclinic with parameter 0=0° appeared as a dark line on the 
screen and was traced on the paper. Both nicols were then rotated 1° 
and the corresponding isoclinic was traced. This process was re- 
peated until the complete pattern of isoclinics shown in figure 2 was 
obtained. 

3. LOCATION OF ISOSTATIC LINES 

When the complete family of isoclinics had been traced, the paper 
was removed and the Q isostatics shown in figure 2 were drawn in by 
the method already outlined and illustrated in figure 1 . 

4. EVALUATION OF STRESSES 

The paper screen was again inserted in the optical system, and the 
Babinet compensator was introduced between the model and the 
analyzing nicol. For each determination the axes of the compen- 
sator were turned so as to be parallel with the lines of principal stress 
at the point considered. The nicol prisms, with their axes still per- 
pendicular to each other, were turned so that the axes made positive 
and negative angles of 45° with the vertical. At every point where 
an isostatic was crossed by an isoclinic and at other suitable points, 
measurements were made of the compensator displacement necessary 
to nullify the retardation produced in the model. From these values, 
which are proportional to the differences of principal stresses, as has 
just been shown, the magnitudes of (P— Q)/T were determined by the 
use of eq 15. 

In all this work, when the load was removed the rubber showed no 
evidence of double refraction even after repeated cycles of stress. This 
fact proved the absence of initial residual stresses, the presence of 
which is one of the commonest faults of most model materials. 

Sufficient data were now available for the evaluation of the principal 
stresses by the use of one or more of the eq 4, 5, 6, and 7. For example, 
in the data which were taken along the z-axis and which are given in 
table 2, the isostatic is itself an isoclinic, and, as already stated, this 
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is the ideal case for the application of eq 5. The logical place to 
begin the integration is at a free boundary, where the Q stress is 
known to be zero. Unfortunately, in the present instance no obser- 
vation was made at the boundary. Therefore, the value —0.055, 
obtained by the theoretical solution for an infinite sheet as given in 
the section below, was taken for Q /T at the first point of observation. 
The magnitude of such an assumed value is likely to be slightly too 
large because of the finite size of the sheet. 

At successive points along the isostatic, four observations were 
made: (1) S Q , the distance along the Q stress line from the point of 
commencement of the integration to the point considered; (2) AS P , 
the distance, measured along the P stress line through the point, 
between the two nearest isoclinics lying on opposite sides of the Q 

i.o 

0.8 



0.6 
0.4 
0.2 

■0.2 



-0.4 



sjj**^ 


O 


P/T 

o 
































THEOR 


E TIOA L 




Uj 




• \sl-tl- *-l 






EX PER 


'MENTAL o « 














1 


1 - v- ^-r* r:sr "" * 










\^ •***' 


Q/T 











2 3 4 

DISTANCE RATIO - r/(L 



Figure 3. — Principal stresses along the x-axis. 

Later independent observations failed to indicate as large a discrepancy between experimental and theoret- 
ical values as is indicated for the values of P/T at r/a= 2.97 and r/a= 3.72. Therefore, the experimental 
curve has not been continued through these points. 

line and approximately equidistant from it, this approximate equality 
being brought about, if necessary, by the proper selection of the point 
to be considered; (3) A0, the amount in radians by which these 
isoclinics differ; and (4) (P—Q)jT from the compensator reading. 

Graphic evaluation of the integral in eq 5 could be made by measure- 
ment of the area under the curve whose coordinates w T ere ( — tjt^ W 

and S Q . These curves were not actually drawn in the present work, 
but the value of the integral was obtained for each interval as the 
product of mid-ordinate by the interval length. 

The addition of values obtained in this manner yielded the results 
given for Q e /T in table 2 and in figure 3. Finally, computations 
were made from the relation P/T=(P—Q)/T+Q/T to obtain the 
results shown for PJT. 

Along the ?/-axis, where again the isostatic is itself an isoclinic, 
there occurs the ideal case for the application of eq 4. The operation 
of obtaining P/T here was carried out in quite similar fashion to 
what has just been described in the application of eq 5. The value of 
Po/T was taken from the theoretical solution as 1.054. The results 
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are shown in table 3 in the column headed P e /T. Values of Q/T ob- 
tained from the relation Q/T=P/T—(P—Q)/T are given in table 3 
in the column headed QejT. Since these values are obtained as the 
differences of two quantities which are much larger than these differ- 
ences, the accuracy would not be expected to be so good as if they 
were obtained more directly. 



Table 2. — Principal stresses along the x-axis 






Distance along the Q isostatic, measured in arbitrary units of such size that the radius of the disk 

(a =2.61 cm) is equal to 30 units. 
Distance along the P isostatic, measured in the same units as &,. Since eq 5 and 6 involve merely 

the ratio of these quantities, the size of the unit is immaterial, 
compensator displacement, 
parameter of the isoclinic, namely, the angle from the direction of the applied tension to the direction 

of the P stress, 
distance from the center of the disk, 
radius of the disk, 
applied tension. 

• experimental value of principal stress ratio Q/T, calculated by eq 5. 
theoretical value of principal stress ratio Q/T, calculated by eq 17 to 22. 
^experimental value of principal stress ratio PIT, calculated from Q«/T and the relation P/T= 

(P-Q)!T+Q/T. 
Pt/T= theoretical value of principal stress ratio PJT, calculated by eq 17 to 22. 



a-- 
T- 

Q.IT- 

QtfT- 
P.l T 



s* 


X 


A0/ASV 


r/a 


Q./T 


QtIT 


P./T 


Pi\T 


Units 


cm 


Radians/ unit 


OO 

4.23 
3.72 
2.97 
2.17 

1.77 
1.58 
1.33 
1.23 
1.16 

1.11 
1.09 
1.06 
1.00 




0.000 
-.055 
-.070 
-.106 
-.179 

-.238 
-.261 
-.297 
-.290 
-.271 

-.246 
-.233 
-.207 
-. 133 




1.000 
.988 
.983 
.969 
.924 

.855 

.780 
.618 
.497 
.379 

.272 

.224 

.141 

-.067 


0.0 
15.2 
37.8 
61.8 

73.8 
79.5 
86.8 
89.9 
92.0 

93.5 
94.0 
95.0 
96.8 






"-0.055 
-.064 
-.096 
-.159 

-.216 

-.250 
-.289 
-.296 
-.293 

-.282 
-.272 
-.241 
*-. 170 




1.509 
1.554 
1.749 

1.779 
1.734 
1.499 
1.384 
1.044 

0.834 
.699 
.534 
.419 


-0. 0013 
-.0018 
-.0035 

-.0056 

-. 0060 

-. 0050 

.0000 

.0050 

.0233 
.0700 
.1050 


0.831 
.826 
.878 

.839 

.778 
.600 
.525 
.326 

.212 

.142 

.076 

»-.04 





• Theoretical value. 
b Extrapolated value. 

For the purpose of obtaining values of the principal stresses around 
the circumference of the disk and for obtaining independent values 
of Q/T on the ?/-axis, integrations were carried out along most of the 
Q isostatics shown in figure 2. Along lines A to G eq 5 was used for 
about the first six points and eq 7 for the remaining points. For 
line G' and all those below it eq 7 alone was employed. For those 
points where eq 5 was used, the method of calculation was identical 
with that performed along the x-axis. 

With eq 7 three observations were made: (1) \[/ Q , the angle through 
which the line of stress must be rotated in a counterclockwise direction 
to bring it upon the isoclinic, the value being obtained most accurately 
by measuring the angle between the isoclinic and the horizontal, and 
subtracting from this angle the parameter of the isoclinic; (2) <j>, the 
parameter of the isoclinic ; (3) (P—Q)/T from the compensator reading. 
The integrations here also were carried out by the multiplication of 
mid-ordinates and interval lengths. In all these integrations, except 
along the x-axis, the value of Q /T was taken as zero. P/T was cal- 
culated in the same manner as was done along the x-axis. Some of 
the results obtained at the last points of each integration are given 
in table 4 in columns headed Q e /T and PJT, and the others in the 
columns headed Q e t/T and P e t/T in table 3. 
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Table 3. — Principal stresses along the y-axis 

S r = distance along the P isostatic, measured in arbitrary units of such size that the radius of the 

disk (a=2.61 cm) is equal to 30 units. 
(Sq= distance along the Q isostatic, measured in the same units as & P . Since eq 5 and 6 involve merely 

the ratio of these quantities, the size of the unit is immaterial. 
X— compensator displacement. 
0=parameter of the isoclinic, namely, the angle from the direction of the applied tension to the direction 

of the P stress. 
r«= distance from the center of the disk. 
a m radius of the disk. 
T= applied tension. 
Q./T= experimental value of principal stress ratio Q/T, calculated from PJT and the relation QJT**PjT— 

(P-Q)/T. 
Q,</T= experimental value of principal stress ratio Q/T, calculated by eq 7 with Q o /T=0, and with separate 

integrations along isostatics R, O, M, L, K, J, and / in figure 2. 
Qt/T= theoretical value of principal stress ratio Q/T, calculated by eq 17 to 22. 
P t /T= experimental value of principal stress ratio P/T, calculated by eq 4. 
P,i/T= experimental value of principal stress ratio P/T, calculated from Qei/T, and the relation P/T= 

(P-Q)/T+Q/T. 
Pt/T= theoretical value of principal stress ratio P/T, calculated by eq 17 to 22. 



Sp 


X 


A0/AS Q 


r/a 


QJT 


Qei/T 


Qt/T 


Pe/T 


Pei/T 


Pi/T 


Units 


cm 


Radians/unit 


CO 

4.98 
3.58 
3.02 
2.90 

2.71 
2.56 
2.12 
1.85 
1.69 

1.54 
1.35 
1.19 
1.095 
1.073 

1.060 
1.033 
1.000 






0.000 
-.005 
-.008 
-.008 
-.007 

-.006 

-.004 

.007 

.028 

.052 

.089 
.180 
.330 
.487 
.531 

.564 
.644 
.733 






1.000 
1.054 
1.098 
1.139 
1.15 

1.17 
1.19 
1.26 
1.32 
1.37 

1.42 
1.48 
1.52 
1.51 
1.51 

1.50 
1.48 
1.47 


0.0 
42.2 
59.0 
62.6 

68.3 

72.8 
85.8 
94.0 
98.8 

103.3 
109.0 
113.8 
116.6 
117.3 

117.7 
118.5 
119.5 










« 1. 054 
1.096 
1.146 
1.16 

1.19 
1.22 
1.31 
1.38 
1.42 

1.47 
1.53 
1.56 
1.55 
1.54 

1.54 

1.52 

M.50 


1.18 
1.24 
1.32 

1.41 
1.37 
1.42 


1.773 
1.880 
1.923 

1.959 
1. 975 
2.046 
2.110 
2.126 

2.169 
2.074 
1.971 
1.808 
1.761 

1.741 
1.618 
1.320 


-0. 0019 
-.0035 
-.0040 

-.0047 
-.0052 
-.0066 
-.0073 
-.0080 

-.0083 

-.0073 

-.0034 

.0083 

.0165 

.0159 
.0270 


0.045 
.031 
.021 

.029 
.045 
.095 
.127 
.164 

.186 
.297 
.389 
.479 
.497 

.502 
.559 
*.62 










0.003 

-.030 

-.012 

.058 

.118 
.140 
.248 















• Theoretical value. 
b Extrapolated value. 

Table 4. — Principal stresses at the circumference of the disk 

B— angle from the direction of applied tension to the line joining the center of the disk with the point 
considered. It is negative in the quadrant chosen. 
Q e /T= experimental value of principal stress ratio Q/T calculated by eq 5 or eq 7 with Q j T=0, and with 

separate integrations along isostatics H, Q', Q, F, E, D, C, B, A, and X in figure 2. 
Qt/T= theoretical value of principal stress ratio Q/T, calculated by eq. 17 to -22. 
P,/T= experimental value of principal stress ratio P/T, calculated from Q e /T and the relation P/T— 

\P-Q)/T+Q/T. 
P%/T— theoretical value of principal stress ratio P/T, calculated by eq 17 to 22. 

0,= experimental value of angle from direction of applied tension to direction of P stress, as observed 

from the isoclinics of figure 2. 
0j= theoretical value of angle from direction of applied tension to direction of P stress, calculated from 
eq 17, 18, 19, and 22. 



e 


Qe/T 


Qt/T 


P,/T 


Px/T 


4>* 


$t 


Degrees 










Degrees 


Degrees 





0.62 


0.734 


1.50 


1.47 








-13 


.45 


.548 


1.46 


1.53 


9 


10 


-23 


.33 


.267 


1.48 


1.57 


7 


8 


-32 


.20 


-.001 


1.50 


1.53 


5 


3 


-43.5 


-.02 


-.287 


1.22 


1.35 


1 


-5 


-53 


-.19 


-.449 


1.07 


1.12 


-4 


-12 


-61.5 


-.29 


-.508 


0.82 


0. 854 


-8 


-19 


-68 


-.30 


-.489 


.46 


.625 


-10 


-24 


-75 


-.30 


-.419 


.36 


.379 


-12 


-30 


-80 


-.25 


-.329 


.27 


.217 


-15 


-33 


-90 


-.17 


-.133 


-.04 


-.066 
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VI. THEORETICAL SOLUTION 

The theory of the problem of the stresses in an infinite plate with a 
heterogeneous circular insertion has been treated by Sezawa and 
Nishimura, 16 who give a general solution in terms of the elastic con- 
stants of the plate and of the insertion. The constants used are X, 
the plane stress constant 17 (not the Lam6 constant of three-dimen- 
sional problems) and n, the rigidity modulus. 18 In the present appli- 
cation the insertion is to be considered as completely rigid. There- 
fore, the elastic constants for it may be taken as equal to each other 
and infinite. This assumption reduces very materially the complexity 
of the general equations given by Sezawa and Nishimura. Further- 
more, the terms involving the elastic constants of the plate may be 
expressed as functions of Poisson's ratio 77 by the following equation 19 
from the general theory of elasticity; 

(16) 



X+2 M 



As previously stated, Poisson's ratio for rubber has been found 20 to 
be 0.500 for elongations in the present range. 

Under these conditions the general equations of Sezawa and Nishi- 
mura reduce to the following: 



fr 1L . la 2 , r n , 12a 2 9 aH _J ,,_. 

00 1L la 2 , r, 9 a*! 1 , 1C . 

fe If, 6a 2 , 9a 4 "| . .. , in . 

T = -2L 1 ~5P+57 4 J Sin2e (19 > 

Here is the angle between the direction of applied stress T and the 
line r from the center of the disk to the point under consideration; a 

is the radius of the disk; rr and rd are, respectively, the radial and 

tangential stresses across a section normal to r, while 60 is the stress 
normal to a section parallel to r. 

The general theory of plane stresses yields equations 21 giving stresses 
referred to any axes rotated by a given angle with respect to selected 
axes. These equations can be applied to obtain the principal stress 

ratios PIT and Q/T in terms of the stress ratios frjT, 60/T, and fd/T. 
The angle of rotation of axes, namely, the angle between r and P, is 
the difference between <£, the angle from the direction of applied stress 

18 K. Sezawa and Q. Nishimura. Stresses in a stretched plate containing a circular insertion. Rep. Aeron. 
Research Inst., Tokyo Imp. Univ. 6, 25 (1931). 

* 7 E. G. Coker and L. N. G. Filon. A Treatise on Photoelasticity, see pago 115 (Cambridge Univ. Press, 
1931). 

18 See footnote 17. 

t» See page 116 of reference 17. 

so Wilfred E. Thibodeau and Archibald T. McPherson. Pkotoelastic properties of soft, vulcanized rubber. 
J. Research NBS 13, 887 (1934) RP751. Reprinted in Rubber Chem. Tech. 8, 183 (1886): William L. Holt 
and Archibald T. McPherson. Change of volume of rubber on stretching: Effects of time, elongation, and 
temperature. J. Research NBS 17, 657 (1936) RP936. Reprinted in Rubber Chem. Tech. 10, 412 (1937). 

» E. G. Coker. Photoelasticity for engineers. Gen. Elec. Rev. 23, 873 (1920). 
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to P, and 0, the angle from the direction of applied stress to r. In 
this manner one finds: 

T 2\_T^T^Tzos2(<t>-0)\ K } 

£_ir^-J?_ _Jr-ee__l 

T~2lT~ t T Tcos2($-d)] ^ l} 



and 



tan 2(0-0)=-^; (22) 

rr—66 



Equations 17 to 22 can then be used to obtain at any point the 
theoretical values of P/T, Q/T, and 4> in terms of the coordinates (r, 0). 
The results are shown in tables 2, 3, and 4 in the columns headed 
P t /T, Qt/T, and <j> tj and in graphic form in figures 3, 4, and 5. 

VII. DISCUSSION OF RESULTS 

As might be expected, both P/T and Q/T can be seen to have their 
largest values in the region of the circumference of the disk near the 
axis in the direction of the applied load. The largest value of Q/T, 
according to both theory and observation, is found at the intersection 
of the circumference and this axis. The maximum of P/T along this 
axis, however, by both theory and experiment, is found displaced 
somewhat from this point in the direction away from the disk. The 
maximum of P/T around the circumference is also displaced somewhat 
from this point. The maximum of (P—Q)/T at the circumference is 
found roughly halfway around the circumference of the quadrant. 

What is not so obviously expected is the vanishing of the P/T and 
Q/T values on the boundary of the disk. Near the axis perpendicular 
to the direction of applied load the P/T values become negative, 
thereby indicating a slight compression. Q/T also has both positive 
and negative values, changing sign on the circumference about midway 
around the quadrant. 

The maximum value of P is about 50 percent greater than T, the 
applied tensile stress, while the value of Q attains a maximum of the 
order of two-thirds of the applied tension. 

Errors in the observed values, where they occur, may be ascribed 
to one or more of the following factors: (a) Errors inherent in the 
method used — (1) the admitted limitations of precision inherent in 
photoelastic determinations; (2) the difficulties of a graphic method 
of attack; and (3) the necessity of an extensive graphic integration to 
obtain some single points; (6) Errors entering into this particular 
investigation — (1) The difficulty of making photoelastic observations 
at the interface between the disk and the soft rubber, with the con- 
sequent necessity of some degree of extrapolation in a region where 
the stress variation is rapid; (2) the use of a single mean thickness 
instead of the actual thickness at each observed point; and (3) the 
uncertainty arising from the fact that the first observed point was 
not at the edge of the model. 

The theoretical treatment involved assumptions contrary to the 
facts of the actual case in the following respects: (1) The sheet was 
assumed to be infinite in extent in both directions; (2) the deformations 
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were assumed to be infinitesimal; and (3) the effect of the clamps at 
the ends was neglected. 

In view of all these factors involving possible causes of discrepanc} r , 
only the first four of which may be regarded as defects of the photo- 
elastic method, agreement between observed and theoretical values 
is considered very satisfactory on the whole. The qualitative charac- 
ters of the variations are the same in all cases. The locations of all 
critical values are identical, with the exception of Q/T in figure 5. 
The quantitative agreement seems as satisfactory as could be ex- 
pected, except for those values of Q/T in figure 4, which were obtained 
as the difference of two nearly equal quantities, and except for Q/T 
in figure 5. In the latter case, not only were the observations difficult, 
but as the circumference of the disk is approached, there is a very 
rapid variation of Q/T, as can be seen in figures 3 and 4. 

VIII. CONCLUSIONS 

This work outlines a method of experimental procedure and demon- 
strates, by a reasonably close agreement with the theoretical treat- 
ment of a problem of complex stresses, that rubber can be employed 
as a model material in photoelastic work. Unlike most model 
materials it remains free from residual stresses, and, although it 
exhibits comparatively large deformations when stressed, it has 
such a large stress-optical coefficient that good results can be obtained 
with very small applied loads, less than one one-hundredth the size 
of those usually employed. Stress problems in rubber itself, for which 
theoretical treatment is often very difficult, can be readily solved in 
this manner. The application of results obtained with rubber models 
to cases involving steel and other engineering materials should be 
made with caution, however, on account of the difference in the 
values of Poisson's ratio. 
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